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SYMBOLS 
D 
Total drag coefficient: = Zev s 
Fro, i temearar seoci i icient 
Gmetcat Tiftecoetiicient: = ors 


Priv COciiwe1eny @t maximum Lift to drag ratio- Co) = 
Lo 


Lift coefficient at zero angle of attack 


7 
Total moment coefficient: CC, = 
7 Zor SLp 





NOMeENesecOcriiCcient av Zero angie of avback 
Drag; force parallel to relative wind.. (1b) 


ae 
Maximum Lift to drag ratio: ae 
Lo 
yi 
Momenv Of inertia about Y axis, (slug-ft.*) 
Lift; force perpendicular to relative wind. (1b) 


Pitching moment about Y axis positive nose up (lb-ft) 
Reference area (ft.*) 
Deceleration time (sec.) 
Velocity (ft/sec. ) 


4 
(ak 
Velocity for minimum drag: Melee: te (ft/secay 


Aircraft weight (lbs.) 


Body axis with origin at c.g. positive direction 
forward | 


Body axis with origin at c.g. positive direction 
out right wing 


Body axis with origin at c.g. positive direction 
aownward 





Ke UONDPEARXE WP ®€C HWP WW FF WW bw & 


Distance from pilots’ eye to CRT display 
Acceleration; gravity units ceo 
Subscript; initial value 

Altitude (ft.) 

Reference length (ft.) 


Aircraft mass (slug) 


= 7. 
Normal load factor: ma 7 
/ 2 
Dynamic pressure: = apr 
Time (sec.) 
ashen 
Velocity ratio: i be 


-) 


Displacement in X direction (ft.) 


©) 


Vertical gust velocity (ft/se 


Angle of attack (rad) 


XY 


Angle of attack disturbance due to gust: Oe 
Elevator deflection (rad) 

Increment 

Flight path angle (rad) 

Nominal glide path angle (rad) 

Angle about which & linearized 

Frequency (rad/sec.) 

Natural frequency (rad/sec.) 

Air density (slug/ft.?) 

Deceleration parameter: r= a 


Pitch attitude angle, positive nose up (rad) 
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SUMMARY 


A fixed based simulator study of the approach and 
tangential landing of a class of manned re-entry gliders 
was conducted to determine the piloting techniques and vee 
hicle parameters which best define the approach and land- 
ing problem. Tests were conducted on.5 configurations 
typical of proposed vehicles with maximum lift to drag 
ratios of 3 to 5. The problem was investigated from final 
approach to touchdown, in still and turbulent air, simulat- 
ing longitudinal equations of motion only. 

A three phase landing consisting of a steady glide 
final approach, a nearly constant load factor flare maneuver, 
and a 15 to 20 second shallow angle, final glide was the 
preferred technique. The vehicle parameters found to be 
most important were maximum lift to drag ratio and speed 
for minimum drag. A landing analysis based on fitted para- 
bolie drag polara and a econatant load faotor Flare waa 


found to be valid for the vehicles tested. 





INTRODUCTION 


The manned space vehicle capable of atmospheric re- 
emery and aerodynamic maneuvering to a horizontal landing 
is presently conceived as a fixed geometry wing-like lift- 
ing body of low aspect ratio and high profile drag. The 
vehicle may have relatively low wing loading and no pro- 
pulsion in the landing phase. 

2ne approden and landing of Such vehicles is (enerdevere 
ized by high approach speeds, steep aporoach angles, rapid 
eeceleration in smear horizontal flignt, and the probabimae 
of high speed and sink rate at touchdown. These character- 
istics pose a piloting task which demands conceéntravion 
end accuracy, while leaving scant margin for error. Pais 
is particularly true if landing facilities and navigation- 


guidance aids similar to those which accommodate modern 


ao, 


PeceperiOrmance ilehver aircrals are svoe oe uscc. 

Approach and landing technidaues which minimize the 
pllot judgment required, avoiding the usual flare or pulle 
up just prior to touchdown and providing reasonable capa-~ 
Bility of Pilot escape with todays ejection equl pied. 
are proposed in REFS. 1, 2, 3, and 5. The technique has 
been used successfully with the X-15 research airplane 


(REF. 4). 





The objective of this research was to study the tane 
gential landing of a class of re-entry gliders using a 
fixed-base simulator with emphasis on determining the 
pilot and vehicle parameters which best define the approach 
and landing problem. 

Simulated landings of vehicles with maximum lift to 
drag ratios of 3 to 5 were conducted in still and turbulent 
air. The landing performance was analyzed, and an analytical 
procedure applicable to generalized vehicles developed. 

The principle pilot and author was a designated Naval 
Aviator with eight years of military flight experience and 


over 2400 hours in propeller type airplanes. 





EQUIPMENT AND PROCEDURE 


1. Simulator. 

The simulator consisted basically of a fixed base seat, 
instrument panel, CRT runway display, and side controller 
(FIG. 1). The computing equipment included an Electronic 
Associates Pace analog computer, a Libratrol 500 digital 
computer, with an analog-digital voltage converter. Data 
recording equipment included a six channel Visicorder 906C 
Strip recorder and the analog X-¥ plotter. A detailed 
description of the simulator and associated equipment is 
contained in Appendix A. 

The pilot was provided instrument display of the ve- 
hicle's altitude, airspeed, pitch attitude, flight path 
angle, and glide path error. The runway simulation con- 
Sisted of a porthole view of the runway boundaries shown 
in standard perspective. The simulated dimensions of the 
runway were 10,000 feet long by 200 feet wide. Airplane 
pitch altitude cues were represented by vertical motion 
and position of the runway display. Visual range and 
height cues could be interpreted from the perspective 
view of the landing area. Only motion in the plane of 
Symmetry was simulated and no lateraledirectional task 


was assigned the pilot. 





Te side controller was given essentially constant 
Memeceversus displacement characteristics with no cenveria. 
forces. Breakout forces and dead zone were negligible. 

No control system dynamics were simulated. 

Atmospheric turbulence was simulated by introducing a 
Becvenoular or uniform probability density of angle of 
ettack (5 ) disturbances into the airplane equations. The 


disturbances could be converted to an equivalent vertical 


pose velocity oy Dp 
7 7 


A vertical gust of 6.5 fps RMS at the speed for minimum 


rag was used. 


2. Vehicle Dynamics. 
The basic longitudinal equations of motion for a rigid 
body with the standard axis system fixed in the aircraft 


are (REF. 6): 


(2) L=WCIOCE— mV EF = o (tice 

(2) D+ Wile + yt F&F = 0 (Drag) 

(3) SZ y = ye) (Moment ) 
Pao heremacle  eQUuauTens cave. 

(4) gS . Vitter Fr = Voor Aar rx Keen & 


(5) = Vad f= Vere oe 
es ve 


A= 4 (K-X pp) Tere hy 





where the terms sin & and cosxY have been linearized 


using a perturbation about 36 , a nominal flight path angle. 


The aerodynamic forces 


=- {(«,d5) 
D= f (a,c) 
N= f (a, 5) 


were expanded in a Taylor Series including second order 
terms where necessary to represent the re-entry vehicle 


aerodynamics. 


XG 2f 
(Dua Sov Ge, ee aoe sf 


G 26 7 JCA 
(8) be Sov 5 26, 7 Sea 4 54 5 a BoB y 5 ery 





Sn 
0) = LV SL Om, + Sore 4 BEF 4 2 , «5 


Aerodynamic forces were assumed independent of Mach 
number for this experiment since the maximum Mach number 
Simulated was approximately .50. Experiments with higher 
Mach numbers should consider compressibility effects. 

Equations 1 through 6 with 7, 8, and 9 were mechanized 
to simulate the vehicle dynamics in the experiment. A 
detailed description of the programming and the computae- 


tion of the runway display are given in Appendix A. 





3. Vehicle Characteristics. 

In order to evaluate the possible effects of re-entry 
vehicle non-linearities, it was decided to fit as accurately 
as possible the wind tunnel data of typical vehicles, B and 
C (FIGS. 2 and 3). Tests were also conducted on a slightly 
simplified vehicle designated By in which the second order 
aerodynamic coefficients were eliminated. Two additional 
hypothetical vehicles (A and D) of the same general class 
were tested. In each test the vehicle was programmed to 
Simulate the desired characteristics in both a range of 
trimmed and untrimmed flight conditions most appropriate 
to the individual glider. 

The programmed characteristics of the five vehicles 


tested, A, B, B C, and D, are shown in TABLE I and FIG. 4. 


ales 
The characteristics for all vehicles tested are in the 
clean configuration. The lowering of landing gear and 


flaps was not simulated in this investigation. 


4. Simulator Flight Task. 

The flight task was designed to simulate a straight 
in approach and landing on a 10,000 foot runway using 
essentially standard aircraft instruments and visual dis- 
play. In preliminary experiments the problem was started 
with the vehicle at 20,000 feet altitude. However, the 
approach phase with the aircraft initially positioned on 
a nominal approach path was found to be trivial and the 


test runs were shortened to save computer time. 


— = - 





Each test run was commenced in an equilibrium glide 
condition, initially positioned on a nominal glide path 
to the runway at an altitude of 5,000 feet. The initial 
velocity and flight path angle were varied over the range 
of practical flight conditions. 

In each approach configuration, the altitude at which 
the flare or pull-up was initiated was varied until the 
landing performance for that configuration was considered 
optimum by the pilot. A series of at least 5 landings 
for each approach condition was then made in which the 
pilot attempted to execute the best possible approach, 
flare, and landing maneuver. 

The glider trajectory on all test runs was recorded 
by plotting V versus & on the X-Y plotter. Five of the 


following variables were normally recorded as a function 


of time: 
a) VF bh k 
b) &£ CG 
C) z bh @ 


Pilot comments were made and recorded on each test run. 
No minimum speed for touchdown was specified. The 
maximum speed for touchdown was specified by an assumed 


maximum vertical speed of 10 fps. Landing performance 


ee Be > ep - 7 
- =» .) ae 
= 
a ls 2 





was judged on controllability and accuracy for the flare 
maneuver, touchdown accuracy, ability to maintain vertical 
speed at touchdown within assumed limits, and overall 
consistency of the landing maneuver. Additional comments 
were made on the pilot's opinions of the short period vee 
hicle dynamic characteristics where applicable. 

Simulated turbulence with an RMS equivalent to 6.5 fps 
for the vehicle average velocity during final glide was 
used in tests of vehicles A, B, and D. Pilot comments on 


«4 


the effects of turbulence were also made. 





RESULTS AND DISCUSSION 


1. Description of the Landing Maneuver. 

An evaluation of the preliminary experiments validated 
the description (REFS. 1 through 5) of the landing maneuver 
as divided into five parts or phases. 

The first part is a pilot guidance phase probably 
beginning after re-entry during which the pilot maneuvers 
to a high key position near or over the landing field aided 
by a radarevoice link and/or inertial navigation systems. 

The second part is usually described as a visual 
straight in or circling pattern in which the glider is 
maneuvered to intercept the final approach leg at a pre- 
scribed altitude and airspeed. During this maneuver the 
pilot might use TACAN equipment as a primary navaid. 

The third and subsequent parts are the subject of 
this report and will hereafter be referred to as PHASES I, 
II, and III (FIG. 5). 

The three phases are described as follows: 

PHASE I - FINAL APPROACH. The final approach is a 
Steady state glide at nearly constant airspeed and flight 
path angle aimed at a fixed, predetermined geographic point 
short of the runway. 

PHASE II « FLARE. The flare is essentially a constant 
load factor pulleup initiated at a predetermined altitude 


during which the sink rate is reduced to near the maximum 
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allowable for touchdown. The maneuver is completed at a 
typically low altitude and small flight path angle. 

PHASE III «= FINAL GLIDE. The aircraft is flown in 
essentially me"g", constant & flight as excess speed 
and altitude are reduced for touchdown. During this phase 
the pilot has rather precise control over his altitude and 
flight path angle but essentially no control over decelera- 
tion. Landing gear and flaps are normally lowered in this 
phase. 

Another phase which is not included in this work is 


the aircraft roll-out or slide-out after runway contact. 


2. Quality of Simulation. 

The preliminary evaluation of the simulator indicated 
that the combination of instrumentation and visual display 
provided sufficient realism to enable a pilot to perform 
and evaluate the final approach and landing maneuver. The 
absence of depth perception in the runway display was an 
early difficulty solved .by locating the flight path angle 
indicator adjacent to the CRT display where it could be 
included as a peripheral supplement to the limited physical 
depth cues. Accurate landings could then be performed with 
reference to the scope for attitude cues, and with airspeed, 
altimeter, and flight path indicator providing precise cone 


trol of sink rate, 


Ca ee = = | ss 





a 


During the flare, the flight path indicator and altimeter 
meovidead the orimary cues for controlling Che Bligny cara 
curvature in such a way as to arrive at the desired h and & 
for final glide. The CRT display and attitude indicator 


provided the necessary attitude cues. 


3. Selection of Touchdown Speed. 

Since the final approach and landing problem can be 
divided into 3 rather distinct phases, the analysis may 
be divided in the same manner by simply matching boundary 
conditions. It is also convenient to start the problem 
mm reverse since it was found that the touchdown conditions 
can be determined independent of the approach, flare, and 
final glide. 

in general the selection of a minimum ground contact 
speed is not dependent on a conventional “stall speed". 
The re-entry vehicle, because it is most probably a lifte- 
ing body, nas no well defined C,; maximum. Instead the 
minimum touchdown speed might be determined by any one of 
the following: 

ae Wing burfet 

b. Directional stability limits 
ec. Ground clearance Wine. 

d. Maximum Landing gear extension 


e. Minimum trim speed 
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For this experiment it was assumed that there were 
no limitations on the minimum speed except the ability of 
the pilot to precisely control the touchdown. Elevator 
control power was provided in excess of that required to 
trim at the touchdown speed. Pitch damping was incorpo- 
rated in all vehicles to insure satisfactory short period 
dynamic characteristics at minimum speed. Since the long 
period vehicle dynamics are of paramount interest in this 
study, difficulties with short period motions were thus 
virtually eliminated, (See Section 10 on Short Period 
Characteristics). Under these conditions the vehicle 
could be flown to a touchdown at any speed within the 
acceptable limits of vertical speed at touchdown assumed 
(10 fps). 

It was found in the experiment that if flight prog- 
ressed to a speed below a value characteristic of each 
vehicle, the loss of kinetic energy was so rapid that an 
excessive sink rate usually resulted. This indicated that 
the pilot was unable to increase the angle of attack rapidly 
enough to maintain one "g" flight (L=W). This characteris-~ 
tic speed was found to be slightly below the speed for 
minimum drag (Vp). 


Assuming a characteristic parabolic drag polar where 


(10) C= Ds 7 KC 
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the C; for maximum lift to drag ratio is 


(11) * 
GO [4 
a Ss 
Then 4 
K 
(12) a os oe 
and 
/ 
(13) Ee. = 
: la D max. £V KCa, 
Letting 
(14) u- 
e 
(15) w= & 


Substituting in 


as . 
(16) B= ZoVIG, # ZV IKE 
yields F 
MAS altH 
(17) = 


Then with the drag equation (2) 


(18) eh 2 itd ZEA “f 
Jt a  ZEVe } “a 
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| fa 
In this equation the variable u represents a linear con- 
tribution of the profile drag to the deceleration while 


1/u* 


yee a nyperbolic contribution of the inducedwdaraz: 
The constant terms are a combination of gravitational 
acceleration and vehicle characteristics. 

For flight speeds less than Vp the induced drag pre= 
dominates and deceleration rapidly approaches infinity 
as shown in FIG. 6. 

m Simple statistical analysis of the success 
Simulator landings are shown in FIG. 7. Based on this 
empirical data a minimum touchdown speed for a reasonable 
probability of success should be approximately .88 Vp in 


maceabsence of other requirements for 4 nignereienciae 


Speed as previously discussed. 


4. The Maximum Touchdown Speed. 
The maximum allowable touchdown speed is generally 
a structural limitation imposed by the vehicle. The ver- 


tical speed at ground contact depends on both V and re 


i x = Vosveexr 


FIG. 8 shows h versus & where the limiting sink rate is 


assumed to be 10 fps. In practice the airplane will be 
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landed with some small but finite & found to be minus 2 
to 2.5 degrees in this experiment. FIG. 8 shows a dise- 
advantage of vehicles with high Vp, hence high landing 
speeds, requiring relatively smaller & at touchdown. 

It was found with all the vehicles flown in this 
experiment that a maximum sink rate of 10 fps i8 a prac- 
tical limitation and imposes no piloting problems unless 
landings are attempted at speeds below .88 Vpe These 
landings usually occurred unintentionally if the pilot 
leveled the vehicle high and subsequently ran out of kinetic 


energy before he could lower the vehicle to the ground. 


5. Final Glide. 

The final glide (PHASE III) was found best character- 
ized by the time T between the completion of PHASE II and 
touchdown. This "time-to-go" to touchdown is the interval 
needed by the pilot to correct the errors in judgment and/or 
execution which occur during flare allowing him to pre- 
cisely lower the rapidly decelerating glider to the runway 
with an acceptable sink rate. 

Experience with the X-15 airplane (REF. 4) and flight 
Simulator experiments with fighter aircraft in high drag 
configurations (REFS. 1, 2, and 3), indicates that a min- 
imum of 15 to 20 seconds should be allowed in the final 
glide. Reference 3 states that 30 seconds is usually the 


maximum required for any condition. 





16 


In the simulator flights, 10 seconds or less was found 
generally unsatisfactory resulting in unacceptable “scatter" 
in landing performance measured in terms of touchdown posi- 
tion, sink rate, and speed dispersion. 10 to 15 seconds 
was generally found to be marginal resulting in some scatter 
particularly in turbulent air. 15 to 20 seconds was found 
to be quite satisfactory for all vehicles. Times greater 
than 20 seconds were seldom experienced and cannot be 
evaluated. 

Since the minimum touchdown speed is determined inde- 
pendent of the trajectory, integration of the drag equation 
over the final glide allows the prediction of the speed 
required at flare completion for any desired deceleration 
time. | 


Equation (18) in the form 


(19) Me face ene tly 
Jr Des 
where | @ o 
de oz 
fe 


expresses the final glide deceleration for parabolic drag 
assumptions. 
Equation (19) integrated is 


wie <= 
| Ars ss (24a) 4% C 


(20) (0,-0,) Uy, 


a Oo 
aa] 
Oe 





1? 


where 


2 & 2 
k= U-yeu +(Frp+ Bleu 


ZY mm 
O,= Jay FAIZ 
A+tVZu 
ie ZAVZ LL 
aA-yzZ tl 
Computation with equation (20) is unhandy and a simple 


graphical solution might be constructed for a particular 
value of E by plotting we backwards from touchdown using 
| a nominal %y of about minus 2 degrees. The function is 
quite linear in the range of interest and sufficient accu- 
racy can be obtained by plotting only a few points. FIG. 9 
is the approximate solution of equation (20) for the ve- 
hicles tested assuming touchdown at .88 Vp. Results of 

the experiment. show very close agreement with this analysis 
tending to show the validity of the parabolic drag polar 
assumption. FIG. 10 shows a comparison of trimmed aero- 
dynamic characteristics for vehicles B and C with their 
parabolic approximations. For these examples, the region 
of interest for flare and landing lies between Cr of 


approximately .15 and .40. This region shows quite rea- 


sonable agreement. 


7 ae > —- 4a —— _ 
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The deceleration characteristics of re-entry vehicles 
and their landing speeds define the problem in the final 
glide where a rather fixed time must be allowed. Further- 
more, the deceleration is approximately constant. Consider 


then the form of the drag equation for final glide 


21 . 
a dé KOA ei get tls COW ST: 
At LE be Ze 


Substituting the velocity ratio at Vp 


(22) 


Hef fre iat: 


Neglecting the gravitational force for simplicity 


(23) | Se. | 
a 


Thus it can be seen that the vehicles with higher lift to 
drag ratios will have superior deceleration characteristics. 
Next consider (23) integrated between flare completion and 


touchdown. 


(2l4) Vy ot t+ Ve 


But if Viyn equals Aieye| Vp 


(25) g 
a = 2y, , pial MUS 
“fg Jr Ble AG, 





This simplified equation shows the effect of vehicie para- 
meters on final glide and landing speed. The final glide 
characteristics are necessarily determined by botn decel- 
eration, a function of E alone, and the speed for minimum 
drag, @ function of botn wing loading and the drag char- 
acteristics. An example of this is shown in the table 
below. Values shown were calculated for T of 15 seconds 


using FIG. 9. 


MeaiCLeE | & K Cp Vp Vmin Vig 


W 
S 





Note: K and Do are based on the idealized drag polar. 


The above table depicts four of the five vehicles tested. 
Venicle By is identical to B except in the computer program- 
ming. Vehicles B and C are representative vehicles of the 
re-entry class while A and Dien oe nermee vehicles. 

It may be noted that vehicles C and D, although they 
nave @ifferent vehicle parameters, show identical landing 
characteristics. This is possible since Ve is defined by 
3 vehicle parameters while = is defined by 2. Compsring 


yenicles B and C shows the effect of wing loading. Vehicle 


B is superior in maximum lift to drag ratio but must land 
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at the higher speed. Comparing B to A shows the same 
tendency but even more graphically. However, vehicle A 
must be provided the greater starting speed and distance 
from landing area. From the simulator study it could not 
be determined what the effect of the higher deceleration 
and greater landing speeds would be on pilot opinion of 
the vehicle. The measure of success seemed to depend 
primarily on the T. If the vehicle could complete the 
flare with sufficient kinetic energy for approximately 

15 seconds or more in final glide, the landing could be 
successfully completed independent of the relative decel- 


eration and speed characteristics. 


6. Flare. 

The flare maneuver (PHASE II) was found to be best 
characterized by the normal load factor used during the 
maneuver. Although it was found that the flare was essene- 
tially a "closed loop" maneuver in which the pilot senses 
his airplane's flight path curvature and altitude, the 
maneuver tends to be accomplished at nearly constant load 
factor. In this experiment, the pilot was not provided 
an accelerometer, yet repeatedly the flare found best liked 
in terms of controllability and accuracy was about 1.2 to 
1.3 "g". A series of typical time histories of vet 


calibrated in "g" units are shown in FIG. ll. 





fad 


Representing “favorable” or near “optimum"™ pilot 
comments, these data show the consistent average of the 
1.2 to 1.3 "g" flare for 4 out of the 5 vehicles over a 
variety of initial velocity and altitude conditions. 
Moreover, the simulator tests showed that the favorable 
or most preferred flare load factor was actually independ- 
ent of altitude and velocity of initiation, the pilot tend- 
ing to adjust the flare initiation altitude to give the 
preferred load factor over a number of trial runs, and 
without reference to an accelerometer or the data of that 
particular series of test runs. Flares using higher "g"” 
tended to be difficult to control resulting most often 
in misjudged level-off altitudes and subsequent poor land- 
ing or even catastrophic over shooting and premature ground 
contact. No time history data is available on vehicle A, 
but the same tendency was evident in simulated flight 
experience. 

Analysis of the flare maneuver began by considering 


equation (18) 


(26) See _ AIp LE batt sa 
ae oe ue 


together with the lift equation (1) in the form 


(27) ar_ Zz 





Ze 


For of = wh 
cosa = SJ 


equation (27) is 


(28) : Se C77) 
SL TBP, A 





Then 
St ats 2ZEansu +H 


To show the effect of the constant load factor assump- 
tion, a series of trajectories were calculated for comparison 
with actual test trajectories in V versus & space as re- 
corded on the X-Y¥ plotter (FIG. 12). It may be noted that 
the n equal to 1.2 step function approximation to the actual 
"closed loop" load factor profile gives excellent results 
in predicting the actual trajectory. 

It would appear that the higher “g" trajectory would 
be superior since the conservation of kinetic energy for 
use in the final glide is paramount. However, the problem 
of precisely controlling the flight path curvature at the 
higher turning rates, lends favor to the less severe flare. 

Equation (29) reveals that the slope of the trajectory 
in V versus 2 space,which is related to the very important 


kinetic energy loss in the flare, is a strong function of E. 
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To illustrate this, a family of trajectories for vehicles 
A, B, and C were constructed using the load factor 1.2 @ 
consistent with the empirical data (FIG. 13). All vehicles 
are started with the steady glide at 450 fps. The final 
glide condition in each case is minus 2 degrees. Equation 
(29) was used and the non-dimensional u converted to the 
more familiar co-ordinate V. 

It may be observed that the kinetic energy lost during 
the flare increases rapidly as E decreases. Since the 
Steady glide speed is limited by compressibility to about 
475 fps, and the flare completion speed is limited by the 
necessity of a finite deceleration time, there is a defie 
nite limit to the capability of a vehicle to perform the 
landing maneuver using this technique. This was found in 
test experience with vehicle A wherein the T final glide 
was found to approach zero and landing performance was 
unsatisfactory. Once in the final glide the lower E char- 
acteristics (which might be typical of a vehicle with gear 
and flaps extended) does not seem to be a serious limitation 
if sufficient initial velocity can be provided. Use of a 
higher load factor in the flare decreases the kinetic energy 
toss, but again the precision required seems to limit this 
technique. 

Another serious vehicle limitation might be the alti- 


tude lost during the flare. If the flare must be initiated 
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at a relatively high altitude, the accuracy with which the 
pilot can perform the maneuver must suffer. In this re- 
Spect the simulator, even with its "ideal" altimeter and 
limited physical cues, showed this to be the case. As the 
flare altitude increased, the pull-up maneuver tended to 
become "open-loop" until much of the altitude was dissipated. 
At some lower altitude a second or modified flare was 
usually evident. In general,the high altitude maneuver 

was more difficult and less precise. 


To calculate the altitude loss in a flare consider 


(30) yy . 
“4 Yuandsr 
ey 
oT 
Sh Vig LIF 
2 af Some or om 
7. 


With equation (28) 
(32) die 
(7-9 SV 


This calculation may be carried out stepwise together with 
the trajectory calculation. However, by assuming 
sin t= Lert Fang. 


SV AV 
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Then 
(33) a 3 
Ah = Mg "Va. Sig 
A, 
3g (net) 


With equation (33) the altitude lost in the flare for the 


three trajectories of FIG. 12 is 


VEHICLE A 4780 ft. 
VEHICLE B 1398 fts 
VEHICLE C 1620 eh. 


These results are consistent with the experiment and accurate 

enough for planning an approach since the pilot will actually 

vary his pull-up considerably depending on his judgment. 
Simulator experience with the flare maneuver through 

more than about 2,000 ft. showed them to be difficult to 

perform precisely. Flare maneuvers at less than about 600 ft. 

also proved to be aceite but these are seldom encountered 


if the 1.2 "g" pull-up is used. 


7. Final Approach. 

The steady glide flight conditions of the final approach 
phase (PHASE I) can be determined by calculating a desired 
trajectory backwards from touchdown until an equilibrium 
glide flight condition is reached for any given vehicle. 

The eee flight path angle, and rate of descent in 


the final approach will, of course, depend on the lift-drag 
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characteristics of the vehicle. However, the speed is 
bounded by compressibility effects and the maximum accepte- 
able rate of descent probably depends on pilot opinion. 

The latter limitation was not a difficulty in the simulator 


and therefore could not be evaluated. 


8. Synthesis of an Approach. 

An example of how an approach could be synthesized 
is given here and shown in FIG. 14. Vehicle C is used for 
this example. 

It was assumed that the touchdown speed was not other- 
wise limited and .88 Vp was used. It was further assumed 
that 15 seconds was desired in final glide, and the lift- 
drag characteristics were constant (no gear and flaps 
lowered in PHASE III). Final glide angle was minus 2 
degrees. A1.2 "g" flare was planned. 


For these conditions 


Qik 
= = xO 
LE f 





> 


From FIG. 8 


Uj = (LIS (YG = = 728 ps) 


y = 


With equation (29) a stepwise trajectory can then be cone 
structed (FIG. 13). The steady glide conditions are found 


to be as follows: 





(= 402 fps 
eae 7 7 
aye 


SAS eS 
I 


With equation (33) 
Ay = ee ae 


= 


this approach was flown on the simulator and the 
actual envelope of 4 trajectories is also shown in FIG. 14. 
The altitude loss was 900 ft. and T varied between avproxi- 
mately 13 and 17 seconds. 

With this very simple analysis the ability of any 
particular vehicle to perform a planned approach can be 
Quickiy evaluated. Simulator experience has shown that 
an actual approach can be very closely approximated by 
this simple method of synthesis. 

the long period dynamics of the re-entry glider land- 
ing can be characterized by elght variables. These are: 

Le Lift V - Velocity 

D - Drag » «- Reference area 

g - Acceleration of gravity P - Air density 

Oo ~- Flight path angle W - Vehicle weight 
This study has shown that for typical vehicles the lift- 
drag characteristics can be expressed in terms of @ para- 


bolic equation for the landing flight regime. If this 
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is done, the landing can be characterized by 5 non-dimensional 


groups: 
7 
“2 
TE py ; 
K - /LW/K 
fan ve Cat / Ca 


is / 


Zaps (eee 


LE le 
The flight path angle & varies from some steady glide 





value given by 
(34) / Ure 
SIMI, = LE Vee 


to a typically small final glide angle of about minus 2 
degrees. 

The normal load factor "n" is pilot selected and 
essentially constant during the flare maneuver. 

The velocity is normalized by Vp which contains pure 
vehicle parameters. 

The maximum lift to drag ratio E is a pure vehicle 
characteristic, but is independent of wing loading allow- 
ing a degree of freedom between the 2 vehicle parameters 
E and Vp. 

The deceleration parameter 7 is a nonedimensional 
time in which the effects of vehicle parameters on a dee 


sired deceleration time in final glide are shown. 
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PiG. 9 u versus T for final glide, and FIG. 15 for 
the final approach and flare show how the 5 non-dimensional 
parameters define generalized trajectory characteristics. 
From these generalized characteristics the effect of the 
vehicle varameters E and Vz can be shown as in FIGS. 16 
and 17. Increasing E gives lower final approach glide 
angles and Sink rates, less kinetic energy and altitude 
loss in flare, and longer deceleration times in final 
gliae for a given initial velocity. Increasing Ve résulvs 
in increased final approach sink rates, greater loss of 
altitude and kinetic energy in the flare, higher landing 
speecs but a sligntly increased deceleration time in final 


ee oe fOr a given initial velocity. 


(A 


imce the vehicle parameters E and Vip are enough ta 
separate the trajectory characteristics of vehicles which 
reasonably fit the assumption of parabolic drag polars, it 
seems to suggest that E and Ve could be used as a upaci. 
Tor rating vehicles in the manner of flying qualities 
parameters. If a pilot did not object to the short period 
characteristics of a re-entry vehicle, the long period 
dynamics would then become the more important and dominate 
his opinion of any vehicle. Simulator tests with Vehicles 
C and D which have equal E and Vp, but varied wing loading 


and drag polar indicated this was true. In this case no 
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Significant difference in vehicle performance was noted. 
However, the experiment was decidedly limited in this area 
and needs verification. It seems quite reasonable that 
increases in E would always BS welcomed, but that an optie 
mum in Vp might be found since the effects of VR are not 
encountered consistently in all 3 phases of the landing 


as they are with E. 


9. Landing Dispersion. 

The average landing dispersion for the 4 vehicles 
tested was £300 feet based on 68 landings in which data 
was available. The maximum deviation from the mean touche 
down position for each vehicle was as follows: 


Over short 


VEHICLE A +1000 -900 
VEHICLE B 

(including B, ) +1000 -1500 
VEHICLE C +900 - 1000 
VEHICLE D +1000 -900 


Because the touchdown and flare completion speeds 
were predictable for a given final approach condition, it 
was found that the approach glide could be "aimed" initially 
to give the vehicle a high probability of touchdown within 
a small landing area (21500 feet from a predicted touchdown 


point). In the presence of turbulence the probability that 
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the airplane might be “bumped” into an early landing in 
the last 2 or 3 seconds of the final glide was found to 
exist. 

Because of the low altitude, and distance from the 
field (approximately 8,000-12,000 feet) in final glide, 
the pilot may not have a clear view of the landing area. 
For this reason, it is possible that approach zone range 
markers could be an important landing aid to improve touche 


down accuracy. 


10. Short Period Dynamics. 

Although the primary emphasis in this study was on 
long period vehicle motions, it was found in preliminary 
Simulated flights that the basic airplane had unsatisfac-~ 
tory short period dynamics in the landing flight range. 
The basic reeentry vehicle has of itself little or no 
rotary damping. Therefore, without some pitch rate feede 
back, the short period characteristic motion did not pro-~- 
vide the necessary tight attitude and flight path control 
necessary for precise landing. Particularly in turbulent 
air, the unaugmented longitudinal dynamics proved unsatise 
factory. Moreover, it seems unlikely that a re-entry 
vehicle configuration could perform its varied mission 
satisfactorily without extensive longitudinal stability 


augmentation. 
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Using a constant speed approximation, the simplified 


Short period motion is given by the characteristic equation 
2 a 
(35) D+ (E— Me Ma) r= (Met Lee) =0 
 ¢] 


from REF. 6. The dimensional stability derivatives and 
Short period characteristics for Vehicles B and C are 
given in TABLE II for speeds of 350 and 230 fps. Due 
to equipment limitations /4Z was invariant with flight 
condition. 

The damped natural frequency, as approximated, showed 
fairly good agreement with test data. It may be noted 
that the short period characteristics for these vehicles 
without the pitch damper would normally be rated in the 
"poor" to "bad" region in the "thumbprint" (REF. 7). It 
was not particularly desired to overdamp the short period 
motion but the computer limitations dictated fixed values 
of /fe. 

The damped short period characteristics were found 
satisfactory for the landing study. The minor differences 
between the short period characteristics of the various 
vehicles was negligible so far as the landing problem is 
concerned. The same was found with the difference in 
control power, gust sensitivity, and ©, .- Future studies 


might consider these variables in more detail. 





CONCLUSIONS AND RECOMMENDATIONS 


Tne following general conclusions are made for the 
range of parameters and assumptions of this report: 
1. The final approach and landing technique for the re- 
entry may be divided into 3 phases as follows: 
PHASE I - FINAL APPROACH. Steady state glide at con- 
stant airspeed and flight path angle aimed at a fixed 
geographic point short of the runway. 
PHASE II + FLARE. A nearly constant load factor pull- 
up during which the sink rate is retarded while the 
altitude and flight path angle are reduced to near 
Bevel flirnt at low altitude. 
PHASE III - FINAL GLIDE. Deceleration to touchdown 
mcsscntially one “p" fiient. 
2. For the landing maneuver, the re-entry vehicle aero- 
dynamics may be approximated by a parabolic drag polar. 
3. Sor the parabolic drag assumptions, a re-entry vehicle 
approach and landing may be synthesized based on 5 gen- 
eralized, non-dimensional parameters; n, &, u, E, and 7. 
Eme approacn is calculated in reverse from the svoucnasa wm 


to the steady state final glide matching boundary conditions 


between phases. 
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+. if the minimum touchdown speed is not otherwise Limited, 
aespeea Of aporoximately .88 V2 is the minimum landing 
Speed giving acceptable final glide performance. 

a+ A minimum deceleration time in final glide of 10 to 15 
seconds is required for marginal landing performance while 
15 to 20 seconds yields satisfactory performance. 

6. An average load factor of 1.2 to 1.3 "g" was pilot 
selected as “near optimum" for the flare maneuver independ- 
ent of altitude loss. For trajectory analysis, a step 
function or load factor with amplitude of 1.2 "g" gives 
excellent results. 

7e An average landing dispersion of =300 feet and @ maxie 
mum of 1500 feet from the planned touchdown position 
might be expected based on 68 landings with 4 vehicles 
initially positioned on an optimum final glide. 

8. <A maximum sink rate at touchdown of 10 fps is @ vrac- 
tical limitation and was seldom exceeded for an optimum 
approach. Touchdown usually will occur at a small but 
Memive fiight path angle of about minus 2 degrees: 

pee ore Short period dynamics of the vaste: venteleu ae 
generally unsatisfactory for precise landing performance: 
With pitch rate damoing incorporated, the long per tae 
dynamics predominates the landing performance evaluation 


with E and V, as the most significant vehicle parameters. 





The following general recommendations are mede as a 


result of this investigaticn: 


1. Further simulator study to determine the extent to which 
Bm and Ve determine a re-entry vehicle's landing performance 

as seen by the pilot; i.e., make wide variations of the wing 
loading and parabolic drag equation parameters while keeping 


mn and Ve constant... 


2. fFurtner investigation to determine if T final glide is 
a function of E and Vas i1.€., study the effect of cecelera- 
tion, landing speed, and distance from touchdown at flare 


completion. 


De Determine if an average flare load factor of 1.2 to 


1.3 "g" is valid over a wide range of pilots and vehicles. 


4. Determine if limits of E and V3 are defined by avproach 


Brae sink rate, glide angle and altitude lostin {tae 


bee investigate the importance of slope of Une Lii ceca 
on airplane touchdowm attitude and pilot visibility. 

Investigate the short period dynamics to find the limit- 
ing pitch damping necessary for precise longitudinal con- 


or |. 


6. Determine the effect of large airspeed and altitude 


offsets from the optimum final glide conditions. 
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PeVGRaMED VEHICLE CHARACTERISTICS 


| QUANTITY UNITS A By C D 


S LH £t.? 160 136 | 136 | 160 | 160 | 
a | lbs 5 364.0 4000 | 4000 |400C 5025 | 
‘45 : Lbs/ rts Oe ey 29.4 | 29.4% 125.0 32.4 
Z| stug-re?| 3698 | 4520 | 4520 [3698 [3698 
Ay sae 21.17 | 21.05 | 21.05 |Z | 
GS 0 +.10-1) $310 | Saez ee 
Ms Pos. 15725. | dncZ2 Moet ae 
Fe, 814 2573 | 573 |.814 |.8214 
OS .0588 | .0377 | 0564 |.0575 | .0402 
CE Q- 094 0 4590 0 
Cos .0227 | -.0306] .228 |-.158 | -.0349 
@.* 1.195 | 1.752 | 1,22. | 2eO ee 
Des 0 2502 0 +40 © 
CES O --0072| =. 0070 mane 0 
gam ~.187 | 246552| -.0612 |=. 1878) Smo 
as ~.268 | ~.1245| -.1146|-.268 |-.268 
Crk O ~.0524| 0 O O 
ae - 860 |- 1260 | 1260 |- 860 |- 860 | 
Gi V V V V V 


tL 
NOTES: Cig = equivalent rotary damping for Mg = -10 





TABLE If 
DIMENSIONAL STABILITY DERIVATIVES 
AND 
Per eet SHORT PERIOD CHARACTERISTICS 


-—— B VEHICLE C 
| PSs 


= Vee eee ee ros 


Derivative ey a a 350 





| ae 15S | 
: 17 ~10 
| oe , 
: WT ee -25 
| Prequency 
and 
Damping 
Ratio 
ae -115 
(Mo=0) 
Wr 509 
ni 2923 
| (Mge-l0) 
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APPENDIX A 


Details of the Simulation Equipment 
A.1 Mechanization of Equations. 
The Longitudinal equations of motion were simulated 
using the Electronics Associates Pace analog computer and 
@ Libratrol 500 digital computer linked through the analog 
digital voltage converter. The three longitudinal equations 


for body axes plus kinematic relations are: 


oe L/ iz, 
(A-1) YVoer= — 5 COLe eye 7, AO « + Sy sf 


ye LE ay. Af / z 
(A~2) ad Be ~ Gere ae 50S 7 ee 


+ Oy F + Cogs AI H+ feed, | 
(A-3) g= Zov'= As Poo + Ogg & + Ong F + Gag “Sf 
(A-4) Vane ing i Ae ee 
(A-5) KX = Veet KB -V sinh AS 
(A~6) j= A— (x-x,) aah, 


(2-7) ee G-i- AY 





Az 


where the following approximations have been made: 


CoOL AS = Sf 
ee 


The schematic of the mechanization of the above eaqua- 
P2OnS Vs sniown in FIG. Al, 

The aerodynamic coefficients given in TABLE I were 
determined by a best fit to the aerodynamic characteristics 
ga the range of both trimmed and untrimmed flight conditions 
suitable to the problem. 

Dynamic pressure wes computed in the Libratrol 500 
using analog inputs of velocity and altituce anc atmos- 
pneric density programmed in the computer. A standard 


exponential atmosvhere of the form 


ef" 
lee 
was used in the computation. 
Experimental data was recorded on the analog's X= 
plotter and the 6 channel Visicorder 906C strip chart 
recorder having a galvanometer natural frequency in excess 


of 120 cycles per second. The X-Y plotter was normally 





employed to pilot trajectory data such as: 
aes ee ff VS 7, ar fe vs KC 


The 906C recorded time histories of the following variables: 


Ve, &, h, Hi, OC, and ©, 


A.2 Cockpit Display. 

The overall view of the cockpit including fixed base 
seat, instrument display, CRI runway display, and side 
controlier is shown in FIG. 1. The detailed sketcn of the 
imstrument panel and runway display is shown in FIG. &d. 

The runway display, computed in the Libratrol 500 
using analog inouts, shows a runway 10,000 feet by 200 feet 
in standard verspective. Airviane pitch attitude is also 
an input to the CRT. 

The & indicator shows airplane flight path angle in 
degrees. The airspeed indicator was calibrated in feet ver 
second. The attitude indicator read @ in degrees. The <y 
meter was calibrated to indicate glide path error in hundrecs 
of feet. The direction of the error was such that the air- 
vlane must fly toward the needle. The altimeter was cali- 


brated to read in thousands of feet to an altitude of 2,000 


eet, then rescaled to hundreds of feet. The roll attvituce 


}-~* 


Qu 


ndicator was not used in this study. All indicators were 


CN 


caled for inputs from the analog computer. 





AL 


A.3 Computation of Runway Display. 
The runway display was computed in the Libratrol 500 


using inputs of h and x from the analog comouter. The 


¢y 


Mle eC, Clee percpeCtUiyemcan be developed considerin 


PiG, A3 
Computing the front and rear of the runway: 
Yan 2 = 4, 
Tat 8, = Oxrh) 


Then 


4 
(<8 ) = + Jeet (f, = Fixad) 


% ha 


(A-9) Sp = + Tax fh, = + 


Computing the slope of the runway edges: 











74 
Sav A — a 
ee 
2 2. 
Le) tof 
7c 
y je va. 
Tnen 
Lo 
(A-10) Yew p= a7 = 
a, 





A5 


Using programmed runway dimensions and eye to screen dis-= 
tance, the Libratrol 500 computed the perspective view 
with equations A-8, A-9, and A=-10. An input of airplane 
attitude @& to the CRT completed the display. 

A.4 Simulation of Atmospheric Turbulence. 

The atmospheric turbulence was basically a uniform 
probability distribution of angle of attack with adjustable 
RMS amplitude and bias. The random number capability of 
the Libratrol 500 digital computer was used to generate a 
square wave signal F(t) with finite pulses, equal spacing, 
random amplitude and a frequency of 12 cycles per second 
as shown in FIG. A4(a). The probability distribution is 


shown in FIG. Au(bd). From REF. 8 the power spectral den- 


Ct ae 


as shown in FIG. Au&(c). 


sity Yt) is 





As shown by schematic FIG. A5(a), the random signal 
F(t) was filtered and used to modulatea 1KC carrier, then 
recorded on a Sony 262 C AM tape recorder. The corner 
frequency of the first order filter was 1 radian per second. 
The random signal from the tape recorder was then 


demodulated , biased for equal positive and negative a and 





summea with Xx to give the desired turbulence level. The 
eee vet wes NOormaily adjusted CO an equivalent RMS 
vertical gust of 6.5 fps at the speed for minimum crag. 
The RC ecnstant of the demodulator was .5 seconds. The 
filtered and biased turbulence vower svectral density ark 


vrobability density are shown in FIG. A6. 
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